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Abstract 

Recent results from the PICASSO dark matter search experiment at SNO- 
LAB are reported. These results were obtained using a subset of 10 detec- 
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tors with a total target mass of 0.72 kg of ^^F and an exposure of 114 kgd. 
The low backgrounds in PICASSO allow recoil energy thresholds as low as 
1.7 keV to be obtained which results in an increased sensitivity to interac- 
tions from Weakly Interacting Massive Particles (WIMPs) with masses below 
10 GeV/c^. No dark matter signal was found. Best exclusion limits in the 
spin dependent sector were obtained for WIMP masses of 20 GeV/c^ with a 
cross section on protons of a^'^ = 0.032 pb (90% C.L.). In the spin indepen- 
dent sector close to the low mass region of 7 GeV/c^ favoured by CoGeNT 
and DAMA/LIBRA, cross sections larger than a^^ = 1.41 x 10""^ pb (90% 
C.L.) are excluded. 

Keywords: dark matter, WIMPs, superheated droplets, SNOLAB 



1. Introduction 

PICASSO searches for WIMP scatteringusing superheated liquid droplets, 
a variant of the bubble chamber technique ]l|,|2|. The abundance of ^^F in the 
target hquid C4F10 gives PICASSO an increased sensitivity to spin dependent 
WIMP interactions since, with the exception of neutralino scattering on free 
protons, ^^F is the most favorable nucleus for direct detection. Measurements 
and shell model calculations of nuclear magnetic moments show the spin 1/2 
of ^^F is carried almost exclusively by its unpaired proton, enhancing the 
spin dependent cross section by nearly an order of magnitude compared to 
other frequently used detector materials 0, 0|- The light target nucleus ^^F 
together with the low recoil detection threshold of 1.7 keV render the exper- 
iment particularly sensitive to low WIMP masses below 15 GeV/c^. This is 
especially interesting following the DAMA/LIBRA and recent CoGeNT and 
CRESST results 0,13, 3 which are suggestive of a low mass WIMP solution 
of order 10 GeV/c^. Therefore this work will explore both the implications 
of the new data for searches in the spin dependent sector, and the sensitivity 
to the low mass region in the spin independent sector. Previous results ob- 
tained with the same apparatus at SNOLAB, but using only two detectors 
with higher intrinsic background and with smaller exposure (14 kgd), were 
presented in jsj. 

2. Detector Principle 

The detector medium in PICASSO is an emulsion containing C4F10 droplets 
of about 200 fim diameter in polymerized water saturated acrylamide. Since 
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C4F10 has a boiling temperature of = —1.7 °C at a pressure of 1.013 bar, 
at ambient pressures and temperatures the droplets can be in a moderately 
metastable superheated state. A heat spike created by the energy deposition 
of a charged particle traversing a liquid droplet triggers a phase transition if 
it occurs within a certain critical length (of order tens of nm) and exceeds a 
certain critical energy (of order keV). Both quantities decrease exponentially 
with increasing temperature and are functions of surface tension, latent heat 
of evaporation and superheat, where the latter is defined as the difference be- 
tween the vapor and external pressures of the liquid. Details of the detector 
principle are explained in jo], |l3] . The phase transition is explosive and each 
bubble nucleation is accompanied by an acoustic signal in the audible and 
ultrasonic frequency range, which is recorded by piezoelectric transducers. 

Since the detector captures phase transitions, it performs as an energy 
threshold device which can be controlled by setting the temperature and/or 
pressure. The relation between the energy threshold Eth{T) and the operating 
temperature in C4F10 has been determined by measurements using mono- 
energetic neutron beams and with a emitters of known energies (all at 1 bar). 
The results of these calibrations are shown in Fig. [T] and allow a precise 
description of the temperature dependence of energy thresholds ranging from 
0.9 keV up to 800 keV. Details of these measurements by PICASSO can be 



found in [ll|,112[. The gap in the recoil energy thresholds between 0.9 keV and 
7.6 keV is due to the absence of prominent resonances in the ^^V(p,n)^^Cr 
reaction cross section used for the calibration of the low energy threshold^. 



3. Response to Different Particles 

Since each temperature at a fixed constant pressure corresponds to a 
defined recoil energy threshold, the spectrum of the particle induced energy 
depositions can be constructed by varying the temperature. A summary is 
shown in Fig. [2J WIMP induced recoil energies of ^^F nuclei are expected 
to be smaller than 100 keV and therefore become detectable above 30 °C 
(at 1 bar). Particles which produce only low ionization densities, such as 
cosmic muons, 7 and (3 radiation, become detectable when they create sub- 



^EfForts are ongoing to add points across the gap by using the smaller resonances in 
the '^^V(p,n)^"'^Cr cross section and especially at 5.1 keV by exploiting a resonance in the 
^^F-neutron cross section. In addition a point at 4.7 keV and 42 °C can be inferred for 
C4F10 from neutron calibrations with C4F8 made by other authors fl3t . 
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Figure 1: Calibration curve for the energy threshold of F recoils as a function of tempera- 
ture obtained from measurements with mono-energetic neutrons; a particle measurements 
are shown as open (blue) circles at 21 °C and 25 °C. 



keV energy clusters within sub-nm sized regions; this is only observable above 
50°C (less than 1 keV). These particles are well separated from strongly 
ionizing neutron or WIMP induced recoils, which allows efficient suppression 
of such backgrounds at the level of 10~^ to 10~^°. 

Alpha-emitters produce a different response. In Fig. |2] the a curve with 
the lower threshold energy (higher threshold temperature) was obtained after 
spiking the inactive detector matrix with ^^^Am such that only a particles 
entering the droplets can induce nucleations. At the threshold which corre- 
sponds to a deposited energy of Erfep= 71 keV, only a. particles with energy 
depositions at the Bragg peak trigger nucleation. The higher a energy thresh- 
old shown in Fig. |2] (full dots) was obtained with ^^^Ra spiked detectors. In 
this case the ^^^Ra daughter ^^^Rn diffuses into the droplets and the ^^°Pb 
nucleus with the highest recoil energy in the decay chain (E^ec = 146 keV) 
defines the threshold. As shown in Fig. [H the ^^^Ra and ^^^Am related 
thresholds (open circles) are found to be in good agreement with the energy 
thresholds obtained in the calibrations with mono-energetic neutrons. In 
both cases, if the energy thresholds decrease below E^ep = 71 keV the liquid 
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Figure 2: Response to different kinds of particles in superlieated C4F10. From left to 
right: 1.75 MeV 7-rays and minimum ionizing particles (dot-dashed); ^^F recoils modeled 
assuming the scattering of a 50 GeV/c^ WIMP (red); poly-energetic neutrons from an 
AcBe source (dotted); a particles at the Bragg peak from ^"^^Am decays (open triangles); 
and ^^°Pb recoil nuclei from ^^^Ra spikes (full dots). 



becomes sensitive to smaller dE/dx on the a track. It is important to note 
in Fig. [2], that the response remains flat from 1 — 120 keV. This has been 
confirmed with numerous detectors with large a background and indicates 
that the detectors are within an uncertainty of less than 3% fully sensitive 
to energy depositions above threshold. A more detailed discussion can be 



found in 12 



Since the detectors are fully sensitive to a particles over the entire range 
of the WIMP sensitivity, a particles are the most important background for 
this kind of dark matter search. However the shapes of the WIMP (essentially 
exponentially falling) and of the a (constant) responses differ substantially, 
such that they can be separated by fitting the two contributions (Sect. [H]). 

4. Experimental Setup 

The present PICASSO installation at SNOLAB accommodates 32 detec- 
tor modules. The detectors are installed in groups of four inside thermally 
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and acoustically insulated chambers, serving as a temperature control unit 
with a precision of ±0.1 °C in the range from 20 °C to 50 °C. The current 
detector generation consist of cylindrical modules of 14 cm diameter and 
40 cm height jsj. The containers are fabricated from acrylic and are closed 
on top by stainless steel lids sealed with polyurethane 0-rings. Each detec- 
tor is filled with 4.5 litres of polymerized emulsion loaded with droplets of 
C4F10. The active part of each detector is topped by mineral oil, which is 
connected to a hydraulic manifold in order to allow periodic pressurizations 
of the detectors to reconvert bubbles back into droplets. 

In the most recent detector generation, the emulsion has glycerine and 
polyethylene glycol as the main ingredients. During fabrication the viscosity 
of the non-polymerized liquid is used to suspend the C4F10 droplets ho- 
mogeneously and uniformly. The volume distribution of droplets peaks at 
diameters of around 200 fim. On average the active mass of a detector used 
in this analysis is 90 g of C4F10 corresponding to 72 g of ^^F. The active 
mass is known with a precision of 1% from weighing during fabrication, but 
additional uncertainties might arise due to losses of C4F10 during polymer- 
ization or by diffusion into the matrix. Therefore the active detector masses 
and sensitivities are verified and monitored by measurements with a cali- 
brated AmBe neutron source. The values quoted in Tabled] for the detectors 
used in this analysis are the averages of the ^^F masses determined during 
fabrication and neutron calibration measurements. The quoted errors are: 
individual detector uncertainties from weighing and calibration during and 
after the fabrication process; and a common systematic uncertainty from cal- 
ibrations with the poly-energetic neutron sources in the lab and underground 
(AcBe/AmBe). A description of the fabrication and purification of this type 



of detector can be found in 14 



Each detector module has nine piezo-electric transducers, mounted at 
three different heights on a flat spot, milled into the outside of the acrylic 
container wall. The transducers are ceramic disks (Ferroperm PZ27) with 
a diameter of 16 mm and 8.7 mm thickness and a pressure sensitivity of 
27 jjN/jihax. The piezoelectric sensors are read out by custom made low- 



noise preamplifiers. Details of the electronic read-out are reported in [15 
The trigger has multiplicity one: triggering of any of the nine channels causes 
all channels to acquire data. The trigger is fully sensitive at temperatures 
above 24 °C (recoil energies smaller than 78 keV). One detector in the set 
up, not loaded with active liquid but fabricated and read out as the others, 
serves as a monitor for non-particle related backgrounds such as mine blasts. 



6 



Detector 


Mass 


Exposure 




9{F) 


kg{F)d 


71 


64.66 ± 2.40 ± 1.94 


16.09 ±0.77 


72 


59.87 ± 1.6U ± 1.8U 


17.69 ± 0.71 


131 


82.79 ±3.11 ±2.8 


10.89 ±0.55 


134 


71.61 ±0.80 ±2.15 


15.94 ±0.51 


137 


81.35 ±2.56 ±2.44 


16.33 ±0.71 


141 


68.70 ±2.88 ±2.06 


13.37 ±0.69 


144 


41.51 ±1.60 ±1.42 


6.18 ±0.31 


145 


69.85 ±2.79 ±2.10 


7.83 ±0.39 


147 


66.26 ±2.63 ±1.99 


6.55 ±0.32 


148 


109.53 ±3.27 ±3.3 


3.43 ±0.15 



Table 1: Summary of the performance parameters of all detectors used in this analysis. 
Active masses are normalized to the mass of ^^F present in a module. The quoted mass 
errors are: individual uncertainties from weighing and neutron calibration during and after 
the fabrication process; and a common systematic uncertainty from calibrations with a 
poly-energetic neutron source (AmBe). The indicated values for exposure cover data taken 
over the entire temperature range from 28 °C to 48 °C. 



electronic spikes, detector cross-talks and ambient noise sources. 

The entire installation is surrounded by a 30.5 cm thick water shield, 
which serves as a neutron moderator and absorber. This shielding is made of 
242 cardboard boxes containing square polyethylene bags filled with water, 
with a filling factor of about 75%. At the location of the experiment, a 
depth of 2070 m, 90% of the fast neutrons above 5 keV are produced by {a, 
n) reactions in the surrounding Norite rock, with the remaining 10% being 
fission neutrons. The fast neutron fiux was measured to be ~ 3000 neutrons 
m~^d~^ 16|. In order to estimate the expected neutron fiux reduction by 



the shielding, Monte-Carlo (MC) simulations have been performed which 
included all structural materials, the geometric filling factor of the water 
boxes and self shielding effects due to the presence of other detectors within 
the shielding. The performance of the MC simulation was checked against 
measurements using ^He counters (SNO NCD-counters) and good agreement 
was found. The simulations predict a reduction of fast neutrons from the 
shielding by a factor of 35, giving an estimated neutron induced count rate at 
the level of 1.1 neutrons kg^^d^^ (kg of ^^F) for operation at 5 keV threshold 
energy. 



7 



5. Acoustic Signatures for Background Reduction 



Apart from the different temperature or threshold energy profiles which 
can be used to discriminate different particle interactions in superheated liq- 
uids (Fig. 12]), the acoustic signals themselves can be exploited for the discrim- 
ination of particle and non-particle sources. Calibrations with neutron test 
beams and fast neutrons from AcBe/AmBe sources show that the associated 
waveforms have characteristic frequency and time dependences: a short rise 
time, reaching a maximum after 20-40 /xs, with slower oscillations following 
for several milliseconds. In addition the amplitude distributions of the high 
frequency content (> 18 kHz) of the particle induced wave forms are concen- 
trated in a well defined peak. These features are used to construct variables 
which allow the discrimination of particle induced events from non-particle 
backgrounds. 

Acoustic energy (EVAR): This parameter measures the acoustic energy 
of an event. Frequencies below 18 kHz were found to carry no relevant in- 
formation and are removed by a Butterworth high pass filter applied to the 
Fourier transformed acoustic signal. The waveform is squared and integrated 
over the signal duration, starting from a fixed pre-trigger time. The resulting 
values are then averaged over all active transducers to reduce solid angle ef- 
fects. The resolution at FWHM is ~ 20 % for temperatures tested, while the 
centre of the distribution increases smoothly with temperature. The param- 
eter EVAR is used to define an acoustic energy threshold to stay sufficiently 
away from non-particle related noise signals. Details concerning the under- 



lying physics processes are described in [12|, [ITI]. Since the expected signals 



from WIMP induced recoils have a similar intensity to neutron recoils this 
discrimination variable is of prime importance for dark matter searches with 
superheated liquids. 

Frequency content (FVAR): Studies of the fast Fourier transforms (FFT) 
of particle induced waveforms have shown that the majority of the signal 
power can be found in the frequency range between 20 and 70 kHz. A 
variable FVAR is constructed by taking the logarithm of the ratio of signal 
power in the intervals from 20-30 kHz and 45-55 kHz. This variable allows 
suppression of fractures or secondary events which have a significant deficit 
in signal power in the low frequency window; these events sometimes follow 
a true particle induced event and are caused by a weakening of the matrix. 
Mine blast events are also efficiently removed by cuts applied on this variable. 
A more detailed discussion of the discrimination variables and the event types 
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they are able to discriminate can be found in [8| and especially Fig. 4 therein. 

Signal rise time (RVAR): This parameter reflects the steepness of growth 
of the signals and measures the energy content within the flrst 25 /xs after 
the signal start time. This variable was introduced in order to suppress a 
class of background events with a characteristic slow rise time, but with an 
acoustic energy and frequency content comparable to particle induced events. 
This background became noticeable in detectors with increasingly smaller 
intrinsic a contamination, especially above 40 °C where this background 
increases nearly exponentially. The most probable cause of these events is a 
cascade of secondary vaporizations in the vicinity of primary particle induced 
events. Since these signals have only slightly reduced contribution at lower 
frequencies, they can only be partially removed by the FVAR variable. A 
scatter plot of the variables RVAR vs. EVAR at 45 °C is shown in Fig. [31 
Particle induced events accumulate in the right upper rectangle, secondary 
background events concentrate at low values of RVAR. 
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Figure 3: The signal energy {EVAR) and rise time {RVAR) related variables allow the 
discrimination of particle induced events from other acoustic signals created by activities 
in the detector matrix. Data taken at 45 °C during calibration runs are shown. Neutron 
induced events cluster in the upper right rectangle with well defined EVAR and background 
events concentrate at low values of RVAR. 



9 



6. Data Collection and Analysis 



The analysis presented here was performed on a group of 10 detectors. 
Seven of these detectors were installed at the end of 2008 and belong to 
the most recent generation of PICASSO detectors with the lowest internal 
background and with sufficient exposure to contribute significantly to the 
analysis. The remaining three detectors belong to the previous generation 
and were continuously taking data from June 2007; of these, two are the 
detectors used in the 2009 analysis [8] . A WIMP run typically lasts 40 hours 
after which the detectors are recompressed for 15 hours at a pressure of 6 
bar in order to reduce bubbles to droplets and to prevent excessive bubble 
growth which could damage the polymer. A total of 264 WIMP runs were 
analysed within this period yielding a total exposure of 114.3 kgd in the 
background and signal regions. Approximately every three months calibra- 
tion data have been taken at several temperatures with a weak AmBe neutron 
source (68.71 ± 0.74 s~^), placed equidistant at 10 ± 2 cm from the centre 



of each detector [18[. These data were used to monitor the stability of the 
detectors and to determine cut efficiencies for the discrimination variables 
EVAR, FVAR and RVAR. The combined data from all detector calibration 
runs covering the analysis period are shown in Fig. HI For a given tempera- 
ture, data from all detectors have been combined in a weighted average and 
compared to MC simulations (red curve in Fig. H]). The observed scatter in 
some of the data points is caused by the uncertainty of the source position 
which introduces an additional systematic uncertainty at the level of 5% at 
each temperature point. These tests monitor the long term stability of the 
detectors. They demonstrate that once the count rates have been normalized 
by grams of C4F10 the entire detector array behaves consistently, as one large 
detector. 

The analysis proceeds in the following order: 

• A list of golden runs is established for each detector. To qualify as a 
golden run: at least six working acoustic readout channels are required; 
the duration of the run must exceed 15 h (1 h for calibrations); and the 
gauge pressure in the detector has to be within 0.1 bar of the ambient 
pressure. 

• A selection to remove event bursts with < 3 seconds (< 0.1 seconds 
for cahbrations) between successive triggers is applied. In these low 
background detectors, the probability for successive events within 3 
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Figure 4: Combined data from all detectors from calibration runs with poly-energetic 
neutrons ( AmBe) . Data were taken in regular intervals spread over the entire data taking 
period. For a given temperature, data from all detectors have been corrected for cut 
efficiencies, combined in a weighted average and are compared to simulations (red). The 
threshold energy scale refers to ^^F recoils; for recoiling ^^C nuclei, energies have to be 
multiplied by 1.47. Uncertainties shown are statistical only; the observed scatter of some 
points is due to the uncertainty in the location of the neutron source which introduces a 
systematic uncertainty at the level of 5% at each temperature point 

seconds is negligibly small, and these events are typically retriggers of 
the same events or events physically induced in the detector by primary 
expansion. 

• An event selection is performed on EVAR. This selection is fixed for 
each temperature by fitting a Gaussian curve to calibration data and by 
interpolating between the calibration temperatures. As a large quan- 
tity of bubbles in the matrix leads to decreasing signal amplitude, for 
calibration runs only the first 200 neutron induced events are selected 
in order to maintain acoustic conditions. Selection values are set to 
give 95% acceptance. 

• The events have to pass a selection on RVAR, chosen to yield 95% 
acceptance on calibration data. 
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Detector 144 


30 °C 45 °C 


Triggers/day 
After 3 sec cut 
After EVAR cut 
After RVAR cut 
After FVAR cut 


23.4 ±0.9 60.5 ±1.2 
15.8 ±0.8 40.0 ±1.0 
2.3 ±0.3 3.2 ±0.3 
2.2 ±0.3 2.4 ±0.3 
2.1 ±0.3 2.2 ±0.2 



Table 2: Effect of the applied cuts on the trigger rate at 30 °C and 45 °C. Detector 144 
is shown as an example. 



• Finally the events have to pass a selection on FVAR^ again to yield 
95% acceptance on calibration data. 

The effects of the applied cuts for two temperatures on the trigger rates 
are illustrated for one of the detectors (144) in Table |2j The cut on the 
acoustic energy variable is the most effective discriminator of non-particle 
related signals. The background increases with increasing temperature and 
the cut on the rise-time variable RVAR becomes more important at higher 
temperatures. 

After correcting for cut acceptances and dead time, the events recorded 
by the detectors at each temperature are normalized with respect to the 
active mass (^^F) and data taking time. The count rates of all detectors 
are fiat in the range from 1.7 to 92 keV (48 °C to 25 °C), similar to that 
observed in the presence of a emitters in the droplets (Fig. [2]). The count 
rates averaged over this plateau range are given in Tableland are indicative 
of the level of a contamination in the individual detectors. The decreasing 
rate as a function of detector number reflects the progress in purification 
during fabrication over time. 

The origin of the a background is still uncertain and under investigation. 
It seems probable that a emission occurs within the droplets. This hypothesis 
is supported by studies of the acoustic energy parameter, which show for most 
detectors at 30 °C indications of two groups of events: one characteristic for 
single nucleation by the recoiling a emitter; and a second group of events 
where the recoiling nucleus and a particle add their contributions to the 



acoustic signal. As discussed in |13], at 30 °C this feature is typical for 



detectors where the a emitters are located inside the droplets. Possible 
scenarios are either a direct contamination of the C4F10 itself or diffusion 
of ^^^Rn from ^^^Ra in the polymer matrix into the droplets. Taking as an 
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Detector 


Rate 




l\/f7TitTi 




cts/kg{F)/d 


pb 


GeV/c^ 


71 


327.6 ± 4.3 ± 21.6 


-15.43 ± 8.71 ± 1.4 


10 


72 


134.2 ± 2.9 ± 0.0 


±10.48 ± 7.82 ± 1.0 


9 


131 


31.5 ±1.6 ±2.3 


-1.80 ±3.38 ±0.31 


9 


134 


209.6 ±3.9 ±12.8 


±4.65 ±9.49 ±0.76 


7 


137 


69.9 ±2.1 ±4.7 


±2.76 ± 5.44 ± 0.48 


10 


141 


25.2 ± 1.4 ± 1.8 


-4.71 ±3.53 ±0.19 


12 


144 


60.8 ±3.3 ±4.3 


±1.69 ±6.48 ±0.54 


9 


145 


31.5 ±2.1 ± 12.3 


-0.78 ±5.24 ±0.42 


12 


147 


20.6 ±1.8 ±1.5 


-0.86 ±3.01 ±0.26 


10 


148 


20.0 ±1.9 ±1.3 


-0.28 ±4.30 ±0.33 


8 



Table 3: Summary of analysis results. The averaged rates are corrected for cut efficiencies 
and the systematic errors reflect uncertainties in the mass determination, the detection 
efficiency and the cut-efficiency errors. Cross section values for WIMP interactions on ^^F 
are quoted for a resolution parameter a — 5 (Sect. [7]) at maximum sensitivity of the fits 
obtained for the WIMP mass given in the corresponding column to the right. The sources 
of systematic uncertainties correspond to those listed in Sect. [71 



example detector 148, with the lowest background rate, a contamination at 
the level of 2 x 10~^^ gU g~^ if the activity is located in the C4F10 and of 
2 X 10~^^ gU g~^ for a contamination originating in the matrix is expected. 

In order to combine for illustrative purposes all detectors in a single plot 
of rate vs. threshold energy, the data of individual detectors are renormal- 
ized by their respective a contamination, so that the data can be combined. 
For this the following procedure is carried out: for each detector the av- 
erage count rate over the entire plateau temperature range is calculated 
(28 °C < T < 48 °C); under the hypothesis of absence of WIMPs, this count 
rate is taken as an approximation of the a background level of the detec- 
tor and is subtracted from individual data points at different temperatures; 
data for each detector and temperature are then combined in a weighted av- 
erage; finally temperatures are converted into threshold energies, by taking 
into account that due to the elevated mine pressure (1.2 bar) the measured 
temperature at the location of the experiment corresponds to temperatures 
at surface where the threshold values were calibrated, reduced by 2 °C. The 
resulting threshold energy spectrum shown in Fig. |5] exhibits several note- 
worthy features: the count rates of all detectors as a function of energy are 
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essentially constant; the sensitivity of the experiment for WIMP induced de- 
viations from the constant background is at the level of a few cts kg~^ d~^ 
(kg ^^F); for modest changes in temperature from 28 °C < T < 48 °C the 
dynamic range in threshold energy sensitivity is large and covers the region 
from 1.7—55 keV; errors are dominated by statistics and reflect the time spent 
at respective temperatures; and in terms of sensitivity to light mass WIMPs 
the experiment could still gain substantially by running at the highest tem- 
peratures. Although the background is subtracted here to better visualize 
the spectrum, a fiat background component is included in the overall fit to 
the spectrum during the WIMP analysis. 




Figure 5: Combined data from all detectors for WIMP runs. For each detector the average 
count rate is calculated over the entire temperature/energy range (28 °C < T < 48 °C) 
and subtracted from the individual data points. Data for each detector are then combined 
at each temperature in a weighted average. The rate expected for a hypothetical WIMP 
with Mw = 7 GeV/c^ and (7^^ = 1.2 x 10"'' pb is shown by the red-dotted curve. 



7. Search for a Dark Matter Signal 

To search for a dark matter signal the measured rates as a function of 
threshold energy are compared to those predicted for interactions of WIMPs 
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in our galactic halo with ^^F nuclei, in the presence of a constant a back- 
ground in the detector. We use the formalism described in which approx- 
imates the recoil energy spectrum as an exponentially falling distribution: 

Rq 2/7-1 \ / C2-£'f 



dEp 



where (Er) = 2MfM^/{Mf + Mw^ivw) is the mean average recoil energy; 
Mp and M^y are the masses of the ^^F nucleus and of the WIMP, respectively; 
(vw) is the average velocity of halo dark matter particles and F^{Er) is a 
nuclear form factor taken as 1 for a light nucleus such as fluorine with small 



momentum transfer [19|; and the constants Ci^2 describe the effect of the 
Earth's velocity, Ve, relative to the halo (ci = 0.75, C2 = 0.56 for Ve = 244 km 
s~^). Ro is the expected total WIMP interaction rate per kg of ^^F per day, 

- ^ (55) (osG^) (safe) 

(2) 

where At = 19 is the atomic mass of the target atom; pw is the mass density 
of WIMPs; and ap is the WIMP interaction cross section on ^^F. Since the 
detector operates as a threshold device the observed rate at a given recoil 
energy threshold Er^^{T) is given by: 

P(ER,ER^^{T))^dER, (3) 

where P{Er, Er^^{T)) describes the effect of a finite resolution at threshold 
and the integral extends to Er^^^, the maximum recoil energy a WIMP can 
transfer at its galactic escape velocity of Vesc = 600 km s~^. The shape of 
the threshold curve is discussed in more detail in [ll|, |l2[. It is determined 
by calibrations with neutron sources and a emitters (Fig. ^ and can be well 
approximated by: 

PiER^ER^T)) = 1 -exp (^a(T) (^1 - ^^jf^)) ■ (4) 

The parameter a{T) describes the steepness of the energy threshold. It is 
related to the intrinsic energy resolution and reflects the statistical nature 
of the energy deposition and its conversion into heat. The larger is a, the 
steeper is the threshold. Measurements with a emitters gave a = 10 ± 1 at 
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146 keV. Alpha particles depositing their energy at the Bragg peak seem to 
produce a less steep threshold with a = 5.8 ± 0.7 at 71 keV. Measurements 
by other authors with 17 keV mono-energetic recoils following the reaction 



'^^Cl(nf/j,p)'^^S are compatible with 1 < a < 5 [l3|, but with large uncer- 
tainties. In this analysis we adopt a principal value of a = 5 and let the 
parameter vary within the interval 2.5 < a < 7.5. 

Since for smaller than 500 GeV/c^ the response curves differ in shape 
from the flat a background of each detector, an upper bound on ap is ob- 
tained for each individual detector by fitting the WIMP response curve and 
the fiat a background. For a given Mw the two parameters of the fit are ap 
and a scale factor describing the constant background. The result for each 
detector is shown in Table El Combined in a weighted average, the maximum 
sensitivity occurs for WIMPs in the mass region around Mw = 10 GeV/c^ 
and with ap = —0.72 ± 1.45 ± 0.12 pb (1 standard deviation; a = 5); this 



null result can be converted into a limit [20| on the cross section for WIMP 
interactions on ^^F of ap = 2.00 pb (90% CL.) for resolution parameter 
a = 5. 

The main systematic uncertainties (1 standard deviation) affecting these 
limits on ap are in order of importance: 

• a 3% common systematic uncertainty in the determination of the active 
mass of the detectors, resulting in a 3% uncertainty in the cross-section 
limit; 

• a 3% uncertainty in the recoil detection efficiency inferred from the 
response of a particles; 

• a 2.5% uncertainty in the EVAR cut acceptance and a 1.5% uncertainty 
due to curve fitting of EVAR, results in a 3% uncertainty in the limit; 

• similarly the event selection results in a 3% uncertainty from RVAR 
and a 3% uncertainty from FVAR\ 

• a 1 °C systematic shift in temperature during test beam calibrations 
would result in an energy scale shift, introducing a 1% uncertainty in 
the cross section limits; 

• atmospheric pressure changes at the level of 3% result in uncertainties 

< 1%; 
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• and the hydrostatic pressure gradient of ±2% with respect to the centre 
of a detector module can be translated into an uncertainty of < 1% in 
the cross section. 

The variation of the energy resolution parameter within the uncertainty range 
a = 5 ±2.5 results in a ±1.5% change in the cross section limit at 10 GeV/c^. 
This uncertainty increases at lower WIMP masses and is shown as a broad- 
ening of the limits into confidence bands (Sects. [H]and[n]). 

8. Limits in the Spin Dependent Sector 

The interaction of dark matter particles with nuclei of ordinary matter 
of electro-weak strength has the general form: 

where is the Fermi constant, and Mw,a are the masses of the WIMP and 
detector nuclei respectively [2l|. Ca is an enhancement factor dependent 
on the type of WIMP interaction and F{q^) is a nuclear form factor which 
becomes important for large mass number, A, and momentum transfer, q. 

Spin dependent interactions (SD) with axial vector couplings involve 
squark and Z exchanges and depend on the spin of the target nucleus with 
an enhancement factor of the form: 

= - KiSp) + an{Sn)f (6) 

TT J 

where the effective proton (neutron) coupling strengths, (Sp^n) are the 

expectation values for the nucleon spins in the target nucleus {{Sp) = 0.44 and 



(Sn) = —0.19 in ^^F) and J is the nuclear spin [2l|, |22|, |23|. Assuming that 
scattering of dark matter on ^^F is dominated by interactions with protons, 
the cross section cr^^ for scattering on protons is related to ap by: 

2 (jSD 



Here jjLp^p are the WIMP-proton (fluorine) reduced masses, Cp^ is the en- 
hancement factor for scattering on the free proton and C^^-^ is the corre- 
sponding quantity for scattering on protons in the ^^F nucleus. is 
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obtained by setting a„ = in Eq. and yields the ratio Cp'^/Cp^) = 
1.285 2S, |25|. With Eq. [7] the fit resuh for ap can be converted into a 
cross section on protons of a^^ = -0.008 ± 0.022 ± 0.002 pb (1 standard de- 
viation; a = 5), yielding a best limit of cx^^ = 0.032 pb (90% C.L.) for WIMP 
masses around 20 GeV/c^. The resulting exclusion curve for the WIMP cross 
section on protons as a function of WIMP mass is shown in Fig. |6] together 
with pubhshed results in the spin dependent sector. The broadening of the 
exclusion curve shows the effect of varying the energy resolution parameter 
a within its uncertainty. 
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Figure 6: Upper limits at 90% C.L. on spin dependent WIMP-proton interactions. PI- 
CASSO limits are shown as full lines. Additional curves are from KIMS [26^, COUPP [H 
and SIMPLE [28,j!|. The DAM A/LIBRA [1, [H allowed regions are also shown (Ught 
grey: with ion channelling). Also shown are the spin dependent search results in both 
soft and hard annihilation channels from SuperK [sO] and AMANDA-II/IceCube [sij; and 
theoretical predictions discussed in [32l |33|. 



''The SIMPLE collaboration has recently claimed very competitive limits in 
larXiv:1106.3014( see, however, arXiv:1106.3559l and [arXiv:1107.1515l 
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9. Limits in the Spin Independent Sector 



Spin independent (SI) or scalar interactions proceed via Higgs and/or 
squark exchanges, with Ca in Eq. |5] given by: 

= l^[ZU + {A~Z)U]\ (8) 

where fn,p are the WIMP couplings to the nucleons. For equal couplings 
to neutrons and protons the cross section is proportional to (coherent 
interaction). Using this assumption the cross section becomes: 

(9) 

with A = 19. The limits o\i ap can be translated into an upper bound on the 
WIMP proton cross section in the spin independent sector, with maximum 
sensitivity at Mw = 20 GeV/c^ and a^^ = 6.1 x 10"^ pb (90% C.L.; a = 5). 
The effect of scattering on ^^C nuclei in the target, including a shift in the 
energy threshold, is estimated to be of order 10% and is included in the 
results. 

At the maximum sensitivity these limits are three orders of magnitude 
less stringent than the best limits reached by XENONIOO and CDMS in the 



SI sector [3J,l35| in the range of 50 GeV/c^ < Mw < 80 GeV/c^. However for 
low mass dark matter particles {Mw < 10 GeV/c^) and heavy target nuclei 
the advantage of coherent scattering in SI interactions is largely lost and 
comparable sensitivity can be obtained with a light mass target nucleus, such 
as ^^F, combined with a low energy detection threshold. This low mass region 
has become especially interesting in view of the DAMA/LIBRA and recent 
CoGeNT results which indicate an annual modulation effect for a WIMP 
with a mass of 7 GeV/c^ and a SI cross section close to 1.2 x 10~^ pb js], [cj. 
In the same mass region this analysis excludes cross sections greater than 
a^^ = 1.41 X 10"^ pb (90% C.L.). The CRESST collaboration has also 
reported the observation of an excess of events with a best fit for a dark 
matter particle with a mass of ~ 13 GeV/c^ and a cross section of 3 x 
10~^ pb [tI. Furthermore, this mass range is similar to that required to 
explain the spectrum of 7 radiation observed by FERMI from the galactic 



centre 36 



A summary of allowed regions and exclusion limits in the low mass region 
is shown in Fig. [71 The broadening of the PICASSO exclusion limit is due to 
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the increasing effect of tlie uncertainty in tlie energy resolution parameter, a, 
in tlie low mass region. The interpretation of the DAMA/LIBRA modulation 
effect shown in Fig. [7] in terms of evidence of interactions of dark matter 
particles with ^^Na nuclei assumes a quenching factor of Qnu = 0.3. It 
is interesting to note that this allowed region appears to be disfavored by 
PICASSO using a target nucleus of an atomic weight very close to that of 
22Na. 




Figure 7: PICASSO limits in the spin independent sector (90% C.L.). Only the re- 
gion of recent interest in the range of low WIMP masses is shown. The allowed re- 
gions of DAMA/LIBRA[55, CoGeNT ^ and CRESST and the exclusion limits by 
XENONIOO ^ and CDMS [sH are shown. The broadening of the PICASSO exclusion 
limit is due to the uncertainty in the energy resolution at low threshold energies. 



10. Summary and Perspectives 

The analysis of 10 detectors in the PICASSO set-up at SNOLAB resulted 
in exclusion limits on spin dependent interactions of dark matter particles 
with protons of dp = 0.032 pb at 90% C.L for a WIMP mass of 20 GeV/c^. 
These limits are more stringent by a factor five than the previous PICASSO 
2009 results and with the normal model for WIMP interactions rule out the 
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ion channelling hypothesis invoked to explain the DAMA/LIBRA modula- 
tion effect. The use of the light target nucleus ^^F, combined with the low 
detection threshold of 1.7 keV for recoil nuclei, renders PICASSO particularly 
sensitive to low mass dark matter particles and gives it also some leverage in 
the low mass region of the spin independent sector. The present stage of the 
experiment is approaching the sensitivity to challenge or confirm the claims 
of seasonal modulations by the DAMA and CoGeNT experiments. 

The main improvements with respect to our previous published results 
are: a reduction in a background by up to a factor eight due to improvements 
in detector purification and fabrication; use of a new discrimination variable 
allowing efficient discrimination of non-particle induced events at low recoil 
energy thresholds; and the extension of the analysis from 2 to 10 detectors. 

In the current 32 detector set up eight additional modules have low enough 
background to be used in the standard analysis described here and will be 
included in the analysis, once their exposure gives them sufficient statistical 
weight. Detector modules with higher background will be gradually replaced 
by cleaner modules depending on progress in detector fabrication and purifi- 
cation. 

The implementation of the event by event a recoil discrimination u sing 
the acoustic signal energy discovered by PICASSO and described in [l2|, [l7 
is proceeding and will allow a substantial increase of sensitivity. In order to 
match the anticipated sensitivity of the next stage of PICASSO, the experi- 
ment has been moved to a new location at SNOLAB, allowing an expansion 
of water shielding with a substantial improvement in neutron suppression. 
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